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Constraints on models for the initial collision geometry 
in ultra relativistic heavy ion collisions 
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Monte Carlo simulations are used to compute the centrality dependence of the collision zone ec- 
centricities (£2,4), for both spherical and deformed ground state nuclei, for different model scenarios. 
Sizable model dependent differences are observed. They indicate that measurements of the 2 nd and 
4 th order Fourier flow coefficients V2.4, expressed as the ratio , can provide robust constraints 
for distinguishing between different theoretical models for the initial-state eccentricity. Such con- 
straints could remove one of the largest impediments to a more precise determination of the specific 
viscosity from precision «2,4 measurements at the Relativistic Heavy Ion Collider (RHIC). 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 



Energetic collisions between heavy ions at the Rela- 
tivistic Heavy Ion Collider (RHIC), produce a strongly 
interacting quark gluon plasma (QGP). In non-central 
collisions, the hydrodynamic-like expansion of this 
plasma |1 J results in the anisotropic flow of particles in 
the plane transverse to the beam direction [],[]. At mid- 
rapidity, the magnitude of this momentum anisotropy is 
characterized by the even order Fourier coefficients; 
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14 where 4> p is the azimuthal angle of an emitted particle, 

15 <&rp is the azimuth of the reaction plane and the brackets 
is denote averaging over particles and events. The elliptic 

17 flow coefficient V2 is observed to dominate over the higher 

18 order coefficients in Au+Au collisions at RHIC (i.e. v n oc 

19 (v 2 )% and v 2 « 1) Q, [HJ. 

20 The magnitudes and trends of V2.4 are known to be 

21 sensitive to the transport properties of the expanding 

22 partonic matter [Q Q [], Q _|. Consequently, there is 

23 considerable current interest in their use for quantitative 

24 extraction of the specific shear viscosity, i.e. the ratio 

25 of shear viscosity r\ to entropy density s of the plasma. 

26 Such extractions are currently being pursued via com- 

27 parisons to viscous relativistic hydrodynamic simulations 

28 J 8) , transport model calculations JlJ, [r| and hybrid 

29 approaches which involve the paramctrization of scaling 

30 violations to ideal hydrodynamic behavior |n], [l^, fl3| . 

31 In all cases, accurate knowledge of the initial eccentricity 

32 £2,4 of the collision zone, is a crucial unresolved prereq- 

33 uisite for quantitative extraction of - . 

34 To date, no direct experimental measurements of £2.4 

35 have been reported. Thus, the necessary theoretical esti- 



36 mates have been obtained by way of the overlap geometry 

37 corresponding to the impact parameter b of the collision, 
3a or the number of participants iV par t in the collision zone. 

39 A robust constraint for iVp art values can be obtained via 

40 measurements of the final hadron multiplicity or trans- 

41 verse energy. However, for a given value of -/V part , the 

42 theoretical models used to estimate £2 give results which 

43 differ by as much as ~ 25% Q, Q - a difference which 

44 leads to an approximate factor of two uncertainty in the 

45 extracted 77/s value |i7]| . Therefore, an experimental con- 

46 straint which facilitates a clear choice between the differ- 

47 ent theoretical models is essential for further progress 

48 toward precise extraction of r\j ' s. 

49 In ideal fluid dynamics, anisotropic flow is directly pro- 

50 portional to the initial eccentricity of the collision zone. 



51 A constant ratio for the flow coefficients 



("2) 



0.5 is also 



52 predicted [21| . It is well established that initial eccen- 

53 tricity fluctuations also influence the magnitude of V2.4 

54 si gnificantly ||, [To], | 2^, i.e. the presence of these fluc- 

55 tuations serve to increase the value of i>2,4- Therefore, 

56 one avenue to search for new experimental constraints, 

57 is to use £2,4 as a proxy for U2.4 and study the model 

58 dependencies of their magnitudes and trends vs. N p3 , Tt . 

59 In this communication we present calculated results of 

60 £2.4 for collisions of near-spherical and deformed isotopes, 

61 for the Glauber [[22[ fl | a nd the factorized Kharzeev- 

62 Levin- Nardi (fKLN) J25| |2(J models, i.e. the two primary 

63 models currently employed for eccentricity estimates. We 

64 find sizable differences, both in magnitude and trend, for 

65 the the ratios , £ \.> obtained from both models. This sug- 

66 gests that systematic comparisons of the measurements 
1 for the Ap art dependence of the ratio j^jj for these iso- 
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FIG. 1. Calculated values of £2,4 vs. iVpart for MC-Glauber 
(open symbols) and MC-KLN (closed symbols) for Au+Au 
collisions (a) and near- spherical 148 Dy and deformed 15S Dy 
as indicated in (b). 



topic systems, can give direct experimental constraints 
for these models. 

Monte Carlo (MC) simulations were used to calculate 
event averaged eccentricities (denoted here as £2,4) within 
the framework of the Glauber (MC-Glauber) and fKLN 
(MC-KLN) models, for near-spherical and deformed nu- 
clei which belong to an isobaric or isotopic series. Here, 
the essential point is that, for such series, a broad range 
of ground state deformations have been observed for rel- 
atively small changes in the the number of protons or 
neutrons j2?], ^8). For each event, the spatial distribu- 
tion of nucleons in the colliding nuclei were generated 
according to the deformed Woods-Saxon function: 
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where Rq and d are the radius and diffuseness parameters 
and ^2,4 are the deformation parameters which charac- 
terizes the density distribution of the nucleus about its 
polarization axis (z 1 ). 

To generate collisions for a given centrality selection, 
the orientation of the polarization axis for each nucleus 
(Q\ , <j)\ and 9 2 , <p2 respectively) was randomly chosen in 
the coordinate frame whose z axis is the beam direc- 
tion. For each collision, the values for 7V part and the 
number of binary collisions N co u were determined within 
the Glauber ansatz Q. The associated £2,4 values were 
then evaluated from the two-dimensional profile of the 
density of sources in the transverse plane p s (r±), using 
modified versions of MC-Glauber Q and MC-KLN || 
respectively. 

For each event, we compute an event shape vector S n 
and the azimuth of the the rotation angle 'F* for n-th 



22 harmonic of the shape profile 29 ; 

23 S nx = S n cos (n** ) = /drx/9 s (rj_)o;(rx)cos(n<?!>),(3) 

24 S ny = 



5 n sin(n^*)= / dr±p s (rx)u(r±) sin(n^), (4) 



>F* = -tan' 

n 



Snx 



(5) 



26 where </> is the azimuthal angle of each source and the 

27 weight w(rj_) = rj_ 2 . The eccentricities were calculated 

28 as: 



£2 



(cos2( 



*5)) £4 - (008 4(0-*;)) 



(6) 



30 where the brackets denote averaging over sources, as well 

31 as events belonging to a particular centrality or impact 

32 parameter range. For the MC-Glauber calculations, an 

33 additional entropy density weight was applied reflecting 

34 the combination of spatial coordinates of participating 

35 nucleons and binary collisions IT 7 ! E3| ; 
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where a = 0.14 was constrained by multiplicity measure- 
ments as a function of N par t for Au+Au collisions [po[ . 

The procedures outlined above (cf. Eqs. || - 0) ensure 
that, in addition to the fluctuations which stem from the 
orientation of the initial "almond-shaped" collision zone 
[relative to the impact parameter], the shape- induced 
fluctuations due to nuclear deformation arc also taken 
into account. Note that e 2 .4 (cf- Eq. ||) correspond 
to i! 2 ,4 measurements in the so-called participant plane 
p2| , |24| . That is, the higher harmonic £4 is evaluated 
relative to the principal axis determined by maximizing 
the quadrupolc moment. This is analogous to the mea- 
surement of V4 with respect to the 2 nd order event-plane 
in actual experiments. One consequence is that the den- 
sity profile is suppressed, as well as the moment for the 
higher harmonic. 

Calculations were performed for a variety of isotopes 
and isobars with a broad range of known /3 2 .4 values. 
Here, we show and discuss only a representative set of 
results for 197 Au (R = 6.38 fm, f3 2 = -0.13, /3 4 = -0.03), 
148 Dy (R = 5.80 fm, /3 2 = 0.00, /3 4 = 0.00) and 158 Dy 
(R = 5.93 fm, (i 2 = 0.26, /3 4 = 0.06) §7], H- For thesc 
calculations we used the value d = 0.53 fm. 

Figure |Tj(a) shows a comparison of £2,4 vs. A^part for 
MC-Glauber (open symbols) and MC-KLN (filled sym- 
bols) for Au+Au collisions. The filled symbols indicate 
larger £2,4 values for MC-KLN over most of the consid- 
ered JVp ar t range. The effect of shape deformation is il- 
lustrated in Fig. 0(b) where a comparison of £2,4 vs. 
Ap ar t [for MC-Glauber] is shown for the two Dy isotopes 
indicated. Both £2 and £4 show a sizable increase for 
the isotope with the largest ground state deformation 
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FIG. 2. Comparison of jjfyz vs. -/Vp art for near-spherical 

148 Dy (filled symbols) and deformed 158 Dy (open symbols) 
collisions. Results are shown for MC-Glauber (a) and MC- 
KLN (b) respectively. 
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FIG. 3. iV P a rt dependence of ^ ( a ), ( b ) and ^ 

(c) for Au+Au collisions (see text). The open and closed 
symbols indicate the results from MC-Glauber and MC-KLN 
respectively. 



3 ( 158 Dy). This reflects the important influence of shape- 

4 driven eccentricity fluctuations in collisions of deformed 

5 nuclei | |4| . The magnitudes and trends of all of these 
e eccentricities are expected to influence the measured val- 
a ues of i>2,4 for these systems. 

9 A priori, the model-driven and shape-driven eccentric- 

10 ity differences shown in Fig. [I], need not be the same 
u for £2 and £4. Therefore, we present the ratio vs - 

12 iVp ar t, for both models in Fig. ||. The ratios obtained for 

13 148 Dy (near-spherical) and 158 Dy (deformed) with MC- 

14 Glauber are compared in Fig. |^(a); the same comparison 

15 is given m Fig. |(b) but for MC-KLN calculations. Fig. 
is ||(a) indicates a significant difference between the ratio 
17 for 148 Dy and 158 Dy over the full range of N paxt 
is considered. This difference stems from additional shape- 

19 driven fluctuations present in in collisions of 158 Dy, but 

20 absent in collisions of 148 Dy. The same comparison for 

21 MC-KLN results, shown in Fig. |](b), points to a smaller 

22 difference for these ratios, as well as a different A par t de- 

23 pendence. We attribute this to the difference in the trans- 

24 verse density distributions employed in MC-Glauber and 

25 MC-KLN. 

26 For a given value of A^ par t, the measured ratio of the 

27 flow coefficients t^W for 158 Dy+ 158 Dy and 148 Dy+ 148 Dy 

28 collisions, are expected to reflect the magnitude and 

29 trend of the ratio , £4 w (note that a constant ratio « 0.5 

30 is predicted for ideal hydrodynamics without the influ- 

31 ence of fluctuations Q|). Fig. || suggests that a rela- 

32 tively clear distinction between fKLN-like and Glauber- 

33 like initial collision geometries could be made via system- 

34 atic studies of / v \ 2 for near-spherical and deformed iso- 

35 topes/isobars. Specifically, a relatively smaller (larger) 



36 difference between the ratios 7^32 f° r each isotope, would 

37 be expected for fKLN (Glauber) initial geometries. Sim- 

38 ilarly the scaling of V24 data from the isotopic or isobaric 

39 pair would be expected only for MC-Glauber or MC-KLN 

40 eccentricities. Note that the influence of a finite viscosity 

41 is expected to be the same for both systems and therefore 

42 would not change these conclusions. 

44 The filled symbols in Figs. || (a) and (b) also suggest 

45 a substantial difference in the ratios predicted by 

46 MC-Glauber and MC-KLN respectively, for collisions be- 

47 tween near-spherical nuclei. This difference is also appar- 

48 ent in Fig. ||(a) where the calculated ratios for Au+Au 

49 (P2 = —0.13, /?4 = —0.03) collisions are shown. The 

50 MC-KLN results (filled circles) indicate a relatively flat 

51 dependence for 40 < A part < 200, which contrasts with 

52 the characteristic decrease, for the same N pal - t range, seen 

53 in the MC-Glauber results. 

54 As discussed earlier, each of these trends is expected 

55 to influence the measured ratios of the flow coefficients 
ty*)2 ■ Therefore, an experimental observation of a rel- 

57 atively flat A^p^t dependence for jj^yz [over the range 
sa 40 < Npm-t < 200] , could be an indication for fKLN-like 

59 collision geometries in Au+Au collisions. Such a trend 

60 has been observed in the preliminary and final data sets 

61 reported in Refs. Q an d is consistent with the 

62 conclusions reached in Rcf. [[~[ [+J that the Ap art and 

63 impact parameter dependence of the eccentricity scaled 

64 flow coefficients — and — favor fKLN-like initial collision 

£2 £4 

65 geometries. 

1 The closed symbols in Figs. ||(b) and ||(a) indicate 

2 a decreasing trend for jf^p for near-spherical nuclei for 



4 



Apart > 200. This decrease can be attributed to the 
fact that, in each event, £4 is computed in the reference 
frame which maximizes the quadrupole shape distribu- 
tion, i.e. the so-called participant frame. In this frame, 
£4 can take on positive or negative event-by-event val- 
ues. Consequently, smaller mean values are obtained, 
especially in the most central collisions. Fig. shows 
that the relatively large ground state deformation for 
158 Dy (open symbols) leads to an increase of [rela- 
tive to that for the spherical 148 Dy isotope] which is espe- 
cially pronounced in the most central collisions. However, 
Fig. |||(a) shows that the modest deformation for the Au 
nuc lei docs not lead to a similarly increasing trend for 
N, 



57 and deformed isotopes (or isobars) are required to exploit 

58 these tests. 
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> 200 as implied by data 



17 The relatively flat N part dependence for jj^yz , over the 
is range 40 < 7V part < 200 in Fig. |(a), suggests fKLN-likc 
19 collision geometries. Consequently, it is interesting to in- 
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